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Abstract 

We propose a 3 + 1 dimensional model of gravity which results in inflation at 
early times, followed by radiation- and matter-dominated epochs and a subsequent 
acceleration at late times. Both the inflation and late time acceleration are nearly de 
Sitter with a large hierarchy between the effective cosmological constants. There is no 
scalar field agent of inflation, and the transition from the inflation to the radiation- 
dominated period is smooth. This model is designed so that it yields, at the cost of 
giving up on Lorentz invariance in the gravitational sector, the Dirac-Born-Infeld type 
conformal scalar theory when the universe is conformally fiat. It, however, resembles 
Einstein's gravity with the Gibbons-Hawking- York boundary term in weakly curved 
space-times. 
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1 Introduction 



Scale invariance is ubiquitous in physics and cosmology is no exception. Observations suggest 
that the universe is nearly scale invariant at early and late times, when it is believed to be 
approximately de Sitter space, with the accelerated expansion being driven by a nearly 
constant vacuum energy. 

The approximate scale invariance at early times can be regarded as evidence for cosmic 
inflation [H El [3l IH [5] (see, for more recent reviews, [6l[71IHll9]). Particularly successful is the 
slow-roll paradigm [21 |3l [7j. Not only does it solve fundamental problems of the standard 
Big-Bang cosmology, such as horizon, flatness, and monopole problems, but it also makes 
important predictions, such as a slight tilt from the scale invariant power spectrum of the 
cosmic microwave background (CMB) [lOj. This paradigm typically involves a scalar fleld, 
called inflaton, which acts as the agent of the near exponential expansion of the universe. 
Such an effective theory framework leaves open the lingering question of the particle physics 
origin of the inflaton. Although fundamental approaches, such as String Theory, have abun- 
dant inflaton candidates, (see, for example, [HI [T21 [131 E]), there is hardly any principle 
for model selection within many possibilities. In the absence of a flnal answer to the fun- 
damental nature of the inflaton, it is worthwhile to explore alternatives to inflaton-driven 
inflation. 

In this paper we propose a new model of inflation. Instead of introducing scalar fleld 
agents, we alter Einstein's gravity in a manner that automatically results in inflation at 
early times. Our aim is not to propose the fundamental theory of quantum gravity but 
rather a simple classical model which might well explain the observational data. Our model 
construction is, however, guided by scale invariance and thus the model is, to some extent, 
top-down rather than bottom-up. Interestingly, not only does our model result in inflation, 
but it can also accommodate radiation- and matter-dominated epochs and a subsequent late 
time acceleration with a smooth transition from inflation to radiation. 

Our model is designed so that it yields the Dirac-Born-Infeld (DBI) type conformal 
scalar theory when the universe is conformally flat and resemples Einstein's gravity with 
the Gibbons-Hawking- York boundary term in weakly curved space-times. An uneasy but 
important feature is that it breaks down Lorentz invariance by introducing a preferred time- 
like vector. This is reminiscent of Einstein-aether theory [T3] ; but the time-like vector is not 
an independent degree of freedom from the metric, in our model. Such breakdown occurs in 
the gravitational sector, where observational constraints are less stringent than in the matter 
sector. Moreover, Lorentz invariance is recovered in the weak curvature limit. Our setup 
is also similar to Hofava-Lifshitz gravity [16] in the sense that the breakdown of Lorentz 
invariance means we treat time as distinct from space. 

In section 2 we introduce the model, give our motivation and explain how it is constructed. 
In section 3 we study the Friedman-Robertson- Walker (FRW) universe which results from 
our model. We close our paper with conclusions and discussions in section 4. 
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2 The Model 



We propose the following alteration of Einstein's gravity defined by the action 

(2.1) 



3X 

4^ 



N 



l + ^(i? + /C)-g 



plus radiation and matter actions S^ad + Sm- We work with units c = h = 1 and thus the 
Planck length is ip = \/Gn- There are two dimensionless constants A and q introduced. As 
we will see later, A and q set the scale of inflation and the cosmological constant (CC) at 
late times, respectively. The scalar quantity /C is defined by 

}C:=-^£n (VhK^ = -2 {K^ + rfd^K) , (2.2) 

where h^y = g^y + 11^11^ and K^i, = \£nh^u, with /i, z/ = 0,1,2,3, are the induced metric 
and extrinsic curvature on the space-like surface Ai with the normalised time-like vector n 
normal to M.. 

To be more concrete, we adopt the Arnowitt-Deser-Misner (ADM) decomposition [17] 

ds^ = -N^dt^ + hij [dx' + N'dt) [dx^ + N^dt) , (2.3) 

where N and A^j, respectively, are the lapse and shift functions, and hij with i,j = 1,2,3 
is the spatial metric which coincides with the induced metric defined above. The normal 
derivative is expressed as 

= (dt - £n) , (2.4) 
with the Lie derivative £n for the vector field Ni. The extrinsic curvature then yields 

K^3 = T^r {h^^ - V.A^. - V.AT,) , (2.5) 



2A^ 

where Vj is the covariant derivative with respect to the spatial metric hij. 



A few remarks are in order: First, the action (2.1 ) is related to the DBI action, as we shall 
see shortly, but it is different from Born-Infeld generalisations of Einstein's gravity discussed 
by Deser and Gibbons [18|. Second, it is noteworthy that the theory defined by the action 
(2.1), though generally covariant, breaks Lorentz invariance by introducing the preferred 
time-like vector n. However, in weakly curved space-times, i.e., Gn{R + JC) ^ A, it reduces 
to the Einstein-Hilbert action with the (space-like) Gibbons- Hawking- York boundary term 

mm. 



The motivation to consider the action (2.1) is based on the following observation: with 
the conformally flat ansatz 

ds^ = fp 0(r, xf [-dT^ + dx^] , (2.6) 
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the Einstein-Hibert action (with a cosmological constant) yields the theory 



EH 



WttG 



N 



d^x^^{R + A) 



3 



d X 



+ 



4-3 



(2.7) 



up to total derivative terms. As advocated by Polyakov [21], the IR triviality of the A0^ 
theory might play a role in explaining the smallness of the CCj^ We would rather give 
significance, however, to the fact that this effective scalar theory is classically conformal. In 
fact there is a unique generalisation of the (classical) conformal scalar theory [22] : 



CFT 



-T, 



D3 



d^xA0^ 



'1 + 



A04 



(2.8) 



where A and q are the two arbitrary constants which cannot be fixed by the conformal 
invariance and will be identified with those in (2.1). For a large value of 0, this reduces to 



the A0^ theory. Note that in type IIB superstring theory on AdS^ x S^, this is the DBI 
action for a D3-brane located in AdS^ with as the AdS radial coordinate in the Poincare 
patch, where the constant T^a denotes the D3-brane tension. In the dual M = A super Yang- 
Mills theory, this theory appears as the effective theory for the modulus field = of 
the adjoint scalars in the Bogomolny-Prasad-Sommerfeld (BPS) saturated case q = 1 [23]. 
Incidentally, this is the same action as the one DBI inflation is based on [21]. Note, however, 
that the scalar field is not an infiaton but the scale factor of the metric in our case. 

We first note that, due to the classical scale invariance, the A0'^ theory has the solution 
= ±l/(-\/Ar). Plugging this into (2.6), one obtains de Sitter space. Likewise, the scale 
invariance of the conformal scalar theory (2.8) guarantees the existence of the solution oc 
1/r. As we will see below, the conformal scalar behaves 
0~ 1/(V^ 



~ l/(vAr), near = 0, and 
oo, in the generic case. This has an interesting cosmological 
interpretation, provided we alter the Einstein-Hilbert action so that it yields the conformal 



g-2)r) 



near 



scalar theory (|2.8|) for the conformally flat universe. Namely, the universe will have two stages 

A(l- 



of de Sitter expansion with the CC E: 



2 

inf 



A£p2 



at early times (inflation) and 



q~'^)ip at late times (present acceleration). Choosing g = 1+e with < e ^ 1, we can realise 
a large hierarchy 0(10^'^) between the two CC's. We shall discuss this apparent fine tuning 
in Sec. 4. Upon the inclusion of radiation and matter, the model results in inflation at early 
times, followed by radiation- and matter-dominated epochs and a subsequent acceleration 
at late times. 



The model (2.1) is designed so that it reduces to the conformal scalar theory (2.8), when 

To see it, note first that 



r. 



the universe is homogenous and conformally flat, namely = 
the derivative term ((90)^/0'* inside the square root in (2.8) is almost the scalar curvature, 
since 



R 



-6 



6- 



(2.9) 



^Assuming an appropriate Wick-rotation 

^Observe, however, that this is not so for the generic inhomogeneous conformally flat ansatz. 
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with the ansatz (2.6) and where the 'arrow' speciahzes to the homogeneous ansatz. Thus, 
we need to convert the second derivative of (f) into a first derivative. Herein, we convert only 
the time derivative, at the cost of giving up on Lorentz invariance. It is then straightforward 



to check that (2.2) yields the desired resullj^ as anticipated from the role played by the 



Gibbons-Hawking- York term. 

Let us conclude this section by mentioning a more generic class of models that obeys the 
property discussed in the last paragraph, which has motivated our proposal]^ Consider the 
theory 



N 



G 

1 + -^{R + K: + aWa^ + |3aaa^' + ...)- q 
oA 



(2.10) 



where = n'^V^n'^ is the acceleration of the unit time-like vector n; a, /3 are arbitrary 
coefficients and the dots represent operators of the same type but manifestly non-covariant, 
such as V*ai or with higher spatial derivatives. Since, using the obvious ADM variables for 



(2.6), we have = (0, -6'^dj(f)/(l)^), then 



VX = -^-^^0, (2.11) 

= ^ ^ ; (2.12) 

the addition of these terms is still consistent with our guiding principle. We observe, however, 
that if /3 7^ 0, the model does not resemble Einstein's gravity with a Gibbons-Hawking- York 
boundary term in weakly curved space-times. Thus we take /3 = 0. Taking a 7^ on the 
other hand, yields a potentially interesting model. Indeed, this term is a total derivative and 
therefore does not spoil the desired form of the model in weakly curved spacetime. Moreover, 
adding this term with an appropriate coefficient can yield a Lorentz invariant scalar field 



model for the ansatz (2.6) (but does not restore space-time Lorentz invariance). We shall not 
pursue this more generic class of models further herein. It would be interesting, nevertheless, 
to see if taking a ^ changes the space of solutions of n-DBI gravity found in [27] . 

Finally, we do not wish to break the covariance of the model explicitly. In other words, the 
manifestly non-covariant terms, such as Voj or with higher spatial derivatives are excluded 
from the outset. 



3 FRW universe 



To analyse the advertised remarkable consequences for cosmology of the model (2.1), we 
start by considering the flat FRW model in the conformal frame, i.e., (2.6) with = 0(r). 



^Indeed /C = -6{(j)'^cj) + ^^V^)- 

^We thank the anonymous referee of our paper for raising the possibihty of this generahzation. 
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Then, as emphasised, our model (2.1) reduces to the conformal scalar theory (2.8) in the 
mini-superspace 

3 



S 



1 - 



A04 



1 



(3.1) 



The classical dynamics of the scalar field can be understood intuitively from the D3-brane 
perspective. The second term with g > corresponds to the repulsive force due to the 
Ramond-Ramond 5-form fiux. When q = 1, the repulsive force balances the attractive force 
in the Anti-de-Sitter (AdS) space coming from the first term. This implies that when q > 1, 
the repulsive force dominates. So the D3-brane initially located at the AdS horizon = is 
repelled to the AdS boundary = oo. 



Let us see it more in detail. Since energy is conserved in the field theory (3.1), the 
equation of motion is equivalent to 




S 



(3.2) 



where S is the energy density of this system. Note that from the viewpoint of the scalar 
field dynamics, S is just an integration constant. From the gravity viewpoint, however, this 
is interpreted as the radiation energy. In other words, we need to add the radiation action 



5'rad to the gravity action (2.1) to consider £ ^ 0. This point can be clarified by taking the 



Einstein gravity limit, A — i- oo and g — > 1 with A(l — q) fixed finite, and rewriting (3.2) in 
terms of the scale factor a{t). This yields the Friedmann equation (3.11) with k = Su = 
and the replacement A(l — g^^) by A(g — 1). It then is clear that £ corresponds to the 
radiation energy. 



It is illustrative to rewrite (3.2) as 



1 



02 + vi(j)) = 



with 



V(0) 



-A0^ 



q + 



(3.3) 



The effective potential V{(f)) and the solution 0(r) for g = 1 + e with < e ^ 1 are shown 
in the left and right panels of Fig. [l| respectively]^ There are three distinct regimes: 



(1) Near de Sitter inflation : \/~X(j) <^ ip^ and — (r — Treheat) ^ 



^ ^reheat 

This shows that the inflation scale is E-„yi = \/X£j}, as mentioned earlier. 
(2) Radiation-domination : 1 ^ \f\(^l£^^'^ ^ (1 — q^'^Y^^'^ and r > Treheat 



(3.4) 



^The readers should not confuse this potential with the inflaton potential. 
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Figure 1: (Left panel) The effective potential V^(0) for q = 1 + 10~^ with X = S = 1. In all 
Figures we take ip = 1. (Right panel) The evolution of the universe, log0 vs. conformal 
time T. There are three distinct stages of evolution; (1) Near de Sitter inflation at early 
times, (2) Radiation-domination, and (3) Near de Sitter acceleration at late times. In the 
right panel, the dotted purple, dot-dashed green, and dashed brown lines are the curves 



(3.4), (3.5), and (3.6), respectively. 



28 (t — Tj-eheat) 

where 

Tfeheat IS approximately the time when y/Xcj) ~ £^ 
with the time of reheating. 

(3) Near de Sitter acceleration : VX(f)/£^^^ ^ (1 — 

1 



(3.5) 

and should roughly be identified 



-1/4 



and Trcheat < ^ < 



^cnd 



T - Tend 

Note that this implies the energy scale of the CC to be 



(3.6) 



(3.7) 



much smaller than the scale of inflation. The CC is 10~^^ GeV (which implies the Hubble 
scale Ex ~ 10~^°Mp) and the inflation scale is believed to be about 10^^ GeV, close to the 
GUT scale. Thus in order to generate this hierarchy, we must choose e = g — 1 ~ 10"^^''. 
This clearly requires a fine-tuning (of the same order as in the traditional description). 



Intriguingly, this shows that our model (2.1) results in the flat FRW universe with infla- 
tion at early times, followed by the radiation-dominated epoclijand a subsequent acceleration 



at late times. Both inflation and late time acceleration are nearly de Sitter with a large hi- 

(9(e^^/^), between the effective cosmological constants. There is no 



1-2 



erarchy, l/y 1 — g" 

scalar field agent of inflation, and the transition from the inflation to radiation-dominated 
period is smooth. Note that the very beginning of the universe at r — )■ — oo is exactly de 
Sitter with only a coordinate singularity. (This coordinate patch covers only half of de Sitter 



^It is not essential that there is only the radiation-dominated epoch at this point. As we will see shortly, 
matter can be easily added and the matter-domination follows the radiation-dominated period. 
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space; coordinates may be introduced that extend beyond this singularity [25]. This means 
that, precisely speaking, our flat universe has no beginning of time; it is past eternal.) 

We also emphasise that in this model gravity becomes approximately Einstein's gravity 
from the radiation-dominated epoch onwards. At this period, the equation of motion ap- 
proximates to the Friedmann equation, as will be elaborated below in (3.10) and (3.11). In 
other words, the curvature Gn{R + IC)/{6X) in (2.1) is much smaller than unity. The devia- 
tion from Einstein's gravity becomes significant only in the inflation period. At early times, 
the curvature G'iv(-R + /C)/(6A) ~ —1 and the square root action is very small. As we will 
discuss later, this may indicate that the model (2.1) predicts a rather large non-Gaussianity 
in the spectrum of the CMB. 

We can easily include matter and the spatial curvature. The FRW universe with spatial 
curvature k (not normalised to unit) is given by 



dr^ 



+ r2(de2 ^sin^ dd^"^ 



{3.t 



where the scale factor a{t) and the observer time t are related to the conformal factor 0(r) 
and the conformal time r by 



dt = dp (f){T)dT and a(t) = I 

The effect of the curvature k is only reflected in the scalar curvature R 



(3.9) 



Matter adds the term Sm(P to the RHS of (3.2); this becomes clear once we express the 
equation of motion in terms of the scale factor a{t), yielding a modified Friedmann equation: 



A 



1 + 



Aa2 



1 + 



Aa3 



(3.10) 



where H is the Hubble parameter a ^da/dt. Note that for sufficiently large a, this expands 
to the usual Friedmann equation 



~ A (1 - 



k£ p '2S £p '2S'^£ p 



where we used q 

8nG%p, 



1. It is now clear that we can identify 



a-" 



o 



and 



2£, 



M<-p 



(3.11) 



(3.12) 



The constant ao is the scale factor at present and pc = SuGn/ (^Hq) is the critical density. 
In Fig. [2] the evolution of the flat universe with matter is plotted. In the matter-dominated 
period, the evolution is approximately given by 



(4) Matter-domination : S/Su, (^m/A)^^^ < < (^m/A(1 - q~^)) ' and 



reheat 
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Figure 2: The evolution of the flat universe, log0 vs. conformal time r, with matter for 
g = 1 + 10"'^ with \ = £ = 1 and £m = 10~^. There are four distinct epochs; (1) Near 
de Sitter inflation at early times, (2) Radiation-dominated, (3) Matter-dominated, and (4) 
Near de Sitter late time acceleration. The dotted purple, dot-dashed green, long-dotted red. 



and dashed brown curves are (3.4), (3.5), (3.13), and (3.6), respectively. 



^1 



M 



a ~ 



(3.13) 



where the numerics suggests tm < Treheat- As we remarked in the last footnote, the inclusion 
of matter simply adds the matter-dominated epoch taking over the radiation-dominated 
period. 



Note that comparing (3.10) and (3.11), it is easy to understand why our model behaves 



better, when the scale factor is very small, than Einstein's gravity. For simplicity, let us focus 
on the flat case /c = 0. In the latter, radiation dominates for sufficiently small a, whereas 
in the former, the CC term responsible for inflation is the leading termj^ This is how the 
initial singularity is evaded in our model. In the case k ^ 0, our model still behaves better 
than Einstein's gravity, but it is not completely singularity free, as we discuss next. 

In open and closed universes, the spatial curvature, if not too large, does not change 
much of the evolution except for the early stage of inflation. The evolution of the universe 
with the spatial curvature is shown in Fig. [3} where a slight deviation from near de Sitter 
inflation is visible in earlier stages of inflation. At early times, the scale factor a <^ 1/ -\/A 



and thus the equation (3.10) approximates to 



A 1 + 



Aa2 



(3.14) 



The corresponding effective potentials V{(f)) are plotted in Fig. |4j Note that this is very 



different from the standard case (3.11 ) where the curvature contribution is negligible at early 



^Of course, since quantum gravity fluctuations presumably become important beyond the inflation scale, 
one should not trust the Friedmann equation for very small a. For argument's sake, we are extrapolating 
the Friedmann equation beyond the inflation scale. 
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Figure 3: The evolution of the universe, log(f) vs. conformal time r, for g = 1 + 10^^ with 
X = S = 1, Sm = 10~^, and the spatial curvature \k\ = 10~^. Positive curvature on the 
left panel and negative curvature on the right panel. At early times, the blue line (actual 
evolution) deviates a little from the dotted purple line (de Sitter). For the positive curvature, 
the deviation is downwards, whereas it is upwards for the negative curvature. 



times, comparing to the radiation and matter. So one might wonder if inflation still solves 
the flatness problem in this model. 

The first term drives near de Sitter inflation with the scale i^mf = y/Xij} and the cur- 
vature contribution always becomes significant for sufficiently small a. The positive and 
negative fc's yield faster and slower evolutions, respectively, relative to the flat universe. 
This explains the downwards and upwards shifts of evolution for positive and negative fc's, 
respectively, in Fig. [3j In order for the inflation to occur at all, the curvature scale must 
be well smaller than the inflation scale, k <^ E^„f = . This is the condition that the 



first term dominates over the second term in (3.14) for an appreciable range of a < 1/v A to 



generate the inflation with sufficient e-foldings. Note that the same condition is necessary 
anyway in the standard case too. So, inflation in this model solves the flatness problem in 
the usual manner. 

A few remarks are in order: (1) in the closed universe, in contrast to the flat universe, the 
curvature singularity develops in the infinite past r — — oo. There the curvature contribution 
dominates and the conformal factor goes like 0(r) ~ £p^exp(v^T) which yields ait) ~ ^Jkt. 
Of course, there is no reason to trust this model at the Planck scale and beyond. So this only 
signifies that the absence of singularity, which happens to be the case in the flat universe, 
is not a universal prediction of our model. (2) In the case of an open universe, formally 
speaking, bubbles nucleate from "nothing" , since = is a meta-stable vacuum, as one can 
see from the right panel of Fig. |4| The bubbles of true vacuum inflate at early times and 
then make a smooth transition to the radiation-dominated phase, evading the graceful exit 
problem. Although it is hard to make sense of "nothing", it is noteworthy that this has a 
certain bearing on the "old inflation" of Alan Guth [T]. 
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Figure 4: The effective potential V{(j)) (blue line) near = for g = 1 + 10^^ with A = ^ = 1, 
£m = 10~\ and the spatial curvature \k\ = 10~^. Positive curvature on the left panel and 
negative curvature on the right panel. The dotted purple line is the effective potential for 
k = 0. These show the early time deviations from the flat universe evolution. 

4 Conclusions and Discussion 

We proposed a new model of inflation by altering Einstein's gravity, inspired by the guiding 
principle of scale invariance. With radiation and matter, this model results in near de Sitter 
inflation, followed by radiation- and matter-dominated epochs and a subsequent late time 
near de Sitter acceleration. The transition from inflation to radiation is smooth, and it 
can achieve a large hierarchy between the two effective CC's. Our model breaks Lorentz 
invariance, in the gravitational sector, by the introduction of a preferred time-like vector. 
In weakly curved space-times, however, it resembles Einstein's gravity with the Gibbons- 
Hawking- York boundary term|^ Intriguingly, the Big-Bang singularity is absent in the flat 
univers^ and the open universe realises "old inflation" without the graceful exit problem. 
In a related matter, it should be noted that our model does not appear to have an initial 
conditions problem. 

Similarly to [TB], since our model is endowed with a foliation by space-like surfaces, a 
scalar graviton may remain as a physical degree of freedom [26] . It might be that this scalar 
mode is, in effect, playing the role of inflaton in our model. The Stiickelberg formalism used 
in [26] might be of help to elucidate this point. 

The potential existence of such scalar graviton raises the question if pathologies, of the 
sort discussed in [2B] in the context of Horava gravity, are present in our model. To address 



this issue, we note that the action (2.1), expressing i? + /C in terms of ADM quantities. 



^Observe that even in this limit the symmetry group is different from the full diffeomorphisms group; in 
that sense, our model resembles but is not exactly identical to Einstein's gravity. 

^Even though the energy density does diverge as a 0. This makes the description of the matter content 
in terms of a perfect fluid unreliable; it is still noteworthy, nevertheless, that a diverging energy density 
does not produce a singular geometry. This is analogous, in spirit, to the fact that a point-like charge in 
Born-Infeld electrodynamics, produces a regular electric field at its location. 
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hij, N^, N, as explicitly done in [27], can be rewritten as: 



S 



3A 



/ 



dtd^xVhN \ 1 + 




(i?3 + KijK^i -K^- 2N-^AN) - q 



(4.1) 



N 



where R3 is the Ricci scalar of hij. In particular, observe the existence of the crucial new 
term with the lapse N~^AN. To understand the importance of this new lapse dependent 
term, we note that in [28], a generic action of type 




is considered. The basic observation in [2E] is that if the "potential" V depends on both 
quantities constructed from hij and on quantities constructed from diN/N, then the short 
distance instabilities of the original Horava model are not necessarily present. An analysis 
(to be presented elsewhere) indeed shows that no evidence for such instabilities is found in 
our model. 

From the radiation-dominated period onwards, our model is no different from the concor- 
dance cosmology. However, since it does differ in the inflation period, it may leave distinct 
signatures on the CMB spectra. Due to its near scale invariance in the inflation period, 
the power spectrum in our model will be nearly scale invariant in consistency with the 
Harrison-Zel'dovich-Peebles spectrum at low angular momentum I < 50. It remains to be 
seen, however, whether it can predict a slight red tilt of the spectral index Ug = 0.968 ±0.012 
[To] . Also, as mentioned earlier, our model may predict a rather large non-Gaussianity [29] . 
Indeed, in the inflation epoch, our model deviates significantly from Einstein's gravity, where 
higher derivative corrections become important. As remarked, the DBI square root action 
is small during inflation, so the enhancement of non-Gaussianity may be expected similarly 
to [21]. 
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